Curvature-induced frustration in the XY model on hyperbolic surfaces 
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We study low-temperature properties of the XY spin model on a negatively curved surface. 
Geometric curvature of the surface gives rise to frustration in local spin configuration, which results 
in the formation of high-energy spin clusters scattered over the system. Asymptotic behavior of the 
spin-glass susceptibility suggests a zero-temperature glass transition, which is attributed to multiple 
optimal configurations of spin clusters due to nonzero surface curvature of the system. It implies 
that a constant ferromagnetic spin interaction on a regular lattice can exhibit glasslike behavior 
without possessing any disorder if the lattice is put on top of a negatively curved space such as a 
hyperbolic surface. 
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The minimum energy principle that a physical system 
obeys may be broken by applying external conditions. 
One such condition, which has profound consequences 
in the ergodicity of the system, is geometric confine- 
ment to a curved surface. A typical example is an en- 
semble of classical electrons confined to a spherical sur- 
face [H- It exhibits so many quasistable configurations 
that a unique ground state is hardly observed within a 
feasible time scale [|| . These long-lived states result from 
the incompatibility between crystalline order and surface 
curvature. In fact, partial disorders such as defects and 
dislocations are necessary in order for local crystalline 
order to propagate through a curved surface, wherein 
various possible configurations of disorders cause many 
frustrated states. Such a curvature- induced frustration 
has also been observed on more general geometries 0, Q ■ 

Another important consequence of nonzero surface cur- 
vature is a breakdown of an orientational order. When 
interacting constituents on a curved surface have orien- 
tational degrees of freedom, they can no longer show per- 
fect orientational order. The loss of perfect orientational 
order is due to the noncommutative property of parallel 
transport of vectors On a curved surface, parallel 
transport of a vector along a closed loop docs not main- 
tain its direction but yields a rotation after the round 
trip. This makes it impossible for all vectors to orient 
the same direction, yielding multiple frustrated states at 
low temperatures even if the system does not possess any 
disorder Q . These facts imply a novel class of orienta- 
tional glasses free from any kind of disorder, which is in 
contrast with ordinary spin glasses dominated by certain 
disorder (see, e.g., Refs. @ and Q for attempts at glass 
formation in absence of intrinsic disorder). Understand- 
ing the nature of such curvature-induced glass transition, 
if it exists, should be crucial from viewpoints of statistical 
physics and soft material sciences. Particularly in the lat- 



ter field, systems with curved or fluctuating geometries 
are accessible to synthesize although the interplay 
between the geometry and thermodynamic properties of 
allowed ordered states still remains to be explored @ . 

There are two well-established approaches to analyze 
the orientational order of a physical system: discrete 
lattice simulations and continuum limit approximations. 
The former approach is preferable to study actual evo- 
lutions of orientational order. On curved surfaces, how- 
ever, one cannot construct regular lattices with congruent 
polygons, in general. Hence, the resulting lattice usually 
involves structural defects as mentioned in the first para- 
graph. These defects may give additional contributions 
to the allowed orientational configuration, thus should 
be removed when we are to extract purely curvature ef- 
fects. This can be achieved by employing a surface having 
constant negative Gaussian curvature called a hyperbolic 
surface [3]. This curved surface enables to construct a 
wide range of regular lattices (called a hyperbolic lattice) 
on it, serving as a platform to address the issue. 

In this Rapid Communication, we consider effects of 
curvature-induced frustration on the orientational order 
in the XY spin model defined on a hyperbolic lattice. 
Monte Carlo (MC) simulations are used to evaluate the 
spin configurations at low temperatures revealing the for- 
mation of high-energy spin clusters distributed in the sys- 
tem. We propose that these configurations are multiply 
degenerate due to curvature-induced frustration, which 
may lead to a zero-temperature glass transition as sup- 
ported by calculating the spin-glass susceptibility. 

The hyperbolic lattice we have used is depicted in 
Fig. UK a) in terms of the Poincare disk representation 11 1. 
The lattice consists of equilateral heptagons in the metric 
of the hyperbolic surface so that wc call it a regular hep- 
tagonal lattice. Seven vertices in the central heptagon 
make the first layer denoted by / = 1, which are sur- 
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FIG. 1: (a) A heptagonal lattice with I = 3 in the Poincare 
disk. The circumference of the disk represents the points at 
infinity, (b) Parallelism between two spins in the Poincare 
disk. Parallel transport of an arrow from ith to jth site along 
the geodesic curve (dotted) rotates the arrow by ipij. 
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rounded in a concentric way by the second (I = 2) and 
third (I = 3) layers. The number of vertices N(l) in- 
creases exponentially with Z, making this lattice infinite 
dimensional. This infinite-dimensional property has been 
found to yield various nontrivial thermodynamic proper- 
ties of physical systems defined on this lattice [12j. The 
relevance to the mean-field behavior observed in small- 
world networks has also been discussed [l3| ■ 

Let us imagine two XY spins % and j confined in the 
surface and a geodesic between them. On the analogy of 
the continuum limit, we suppose that i's phase is experi- 
enced by j after parallel transported to j's position along 
the geodesic. The Hamiltonian is then given by 



H 



(1) 



where J is the coupling constant, the sum is over all 
nearest pairs in the system, and (pi represents the spin 
variable at the ith site. Nonzero surface curvature man- 
ifests in the presence of the additional angle xpij. The 
angle ipij describes the amount that i's phase acquires 
after the parallel transport. In the language of differen- 
tial geometry, this additional angle stems from the affinc 
connection [J] between two neighboring sites along the 
geodesic line. Figure [Hb) illustrates how to actually 
determine ipij: we obtain positions of i and j on the 
Poincare disk by hyperbolic tessellation [lj| • If they lie 
on a line through the origin of the Poincare disk, the 
geodesic is represented by a straight line yielding ipij = 0. 
Otherwise, the geodesic appears as an arc of a circle that 
meets the circumference of the Poincare disk at a right 
angle. The positions of i and j together with the origin 
of this circle make an angle, which determines ipij. 

The spins at sites i and j are parallel when fa — <pj — 
ipij = 0. (Note that ipij = —ipji to keep the spin-spin 
interaction symmetric.) Accordingly, parallel transport 
of the ith spin along seven edges of a heptagon alters the 
value of (pi causing frustration in local spin configuration. 
It should be emphasized that the local frustration in our 



FIG. 2: (Color online) Snapshots of spin configurations with 
I = 7, where a darker line represents a higher-energy bond, 
(a) A high-temperature regime with T — 1.00 and (b) a closer 
look at the surface, (c) A low-temperature regime with T w 
0.06 and (d) its magnified image, where clustering structures 
are clearly visible. The bond energy and temperature T are 
in units of J and J/fcs, respectively. 



system has no relation with any quenched disorder and 
originates purely from the intrinsic geometry of the sur- 
face. Such a disorder-free frustration is in contrast with 
disorder-driven frustrations in random gauge XY mod- 
els [13] ■ It is also notable to see how this differs from such 
frustrations as observed in the triangular antiferromag- 
netic Ising ( AFI) model Q : while the hexagonal symme- 
try in AFI allows the presence of loose spins as well as the 
fully-ordered configuration, none of these can be found in 
our system. We point out that the sum of ipij around each 
heptagon is conserved as ^pipij = 27r/ with / = —1/6 
; here J2p represents the summation around a plaquette 
in a counterclockwise direction, and / characterizes the 
strength of frustration. The quantitative invariance of / 
for all constituent heptagons is analogous to that in the 
uniformly frustrated XY model on a flat plane 15| . The 
latter model describes a superconducting film penetrated 
by external magnetic flux [l6j |. and phase ordering in it 
is still under active investigations [17] . It is known that 
the uniformly frustrated XY model exhibits glassy be- 
havior on a square lattice when / is irrational 15, [l8| . 
We should notice, nevertheless, that a direct analogy be- 
tween the two systems is hindered by topological differ- 
ences between planar and hyperbolic lattices. In fact, the 
hyperbolic lattices show glassy behavior for a rational / 
as demonstrated below. 

To equilibrate the system described by Eq. ([1]), we 
employ MC simulations incorporated with the parallel 
tempering (PT) method [l{|. It is a numerical tech- 
nique devised to take averages efficiently over the state 
space. In this method, MC calculations are carried out 
simultaneously at different temperatures, which signifi- 
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cantly saves computational cost for equilibrating large- 
scale frustrated systems. Suppose that we have spin con- 
figurations {4>i}r k and {4>i}T k+1 running at temperatures 
T = T k and T fe+1 , respectively. Without loss of gen- 
erality, we may set Tj~ > Tfc+i- Performing the stan- 
dard Metropolis algorithm on the configurations, we reg- 
ularly check their energies E[{<l>i}T k ] an d E [{0i}T fc+1 ] 
and exchange these configurations with a probability of 
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That is, 



when {4>i}T k+1 is trapped in a local energy minimum, 
PT makes it probe a wider region of the state space by 
passing it to a higher temperature, T k . If it finds a state 
with a sufficiently low energy, the vicinity can be checked 
in more detail by lowering its temperature again. 

Figure [2] plots the resulting spin configurations at high 
and low temperatures. In the former two figures, there 
appears no regularity in the distribution of high-energy 
bonds (denoted by bonds with darker tone) . On the con- 
trary, at low T, high-energy bonds are merged into local- 
ized regions forming clusters. It implies that the forma- 
tion of such high-energy clusters is energetically favor- 
able: even though frustration cannot be removed, they 
may be shifted and confined into small regions, lower- 
ing the energy in other regions instead. Furthermore, we 
have observed that the locations of such high-energy clus- 
ters remain floating by thermal fluctuations. Thereby the 
distribution of high- and low-energy bonds is not unique 
making the ground state degenerate. Indeed, we find that 
each high-energy cluster in Fig. [5] appears as a topologi- 
cal defect, i.e., a negative vortex to reduce the free-energy 
cost induced by the curvature [2(| . It is noteworthy that 
those plots show a similarity to the results in Ref. 2 if . 



which reports that defects of a glass-forming liquid in a 
negatively curved space are concentrated in local regions 
at low T. The origin of negative vortices differs inherently 
from thermal dissociation of a vortex-antivortex pair on 
a hyperbolic surface whose properties were considered in 
Ref. [2^. Still, we note that if the interaction between 
vortices becomes short ranged due to curvature [22| , this 
effect may be relevant in exhibiting glassy features. 

The glassy behavior is quantified by the relaxation 
time to f° r t ne system to evolve from the configurations 
obtained in PT. To this aim, we calculate autocorrelation 
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with the standard Metropolis algorithm for t = 10 4 MC 
time steps. We extract to by fitting C(i) with a stretched- 
exponential function, C{t) oc exp \-(t/t ) b } charac- 
terizing very slow relaxation in glasses. Figures [3ja) and 
[31(b) display the measurements for I = 7 and fitting re- 
sults. The relaxation behavior is described quite well by 
the simple Arrhcnius form, to oc exp(A/T), in which A 
indicates the activation energy. It is noteworthy that the 
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FIG. 3: (Color online) (a) Autocorrelations of spin configu- 
rations, given in Eq. ([2J, together with fitting results with 
a stretched-exponential form. The results are obtained for 
T = 0.250,0.167,0.125,0.100,0.083, and 0.067, from bottom 
to top. (b) Relaxation time to as a function of inverse tem- 
perature from the stretched-exponential fit. Inset: & as a 
function of 1/T from the same fit. It approaches unity at 
high T leading to a simple exponential form. 



persistent linearity of log to with respect to 1/T shown 
in Fig. [3Jb) (i.e., T-indcpcndcnt activation energy A) is 
commonly observed in supercooled liquids such as u-Si02 
and t>-Ge02, in which spatiotemporal fluctuation of local 
configuration leads to critical slowing down 24 1 . 

The presence of a glass transition as well as transition 
temperature can be probed by observing the divergence 
of the spin-glass susceptibility designated by xsg- In 
performing PT to calculate xsg, we prepare two replicas 
fi and v for each T and N . Then we obtain the spin-glass 
susceptibility \sg = NS(T;t' )/T, where 
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with the measurement time t after some equilibration 
time t' [2^ . Figure BJa) gives the T dependence of 
Xsg(T,N) with varying TV's. Our simulation code can 
properly handle the numerical precision up to I = 7, 
and the cases for I < 5 are excluded due to finite-size 
effects. We see that xsg begins to increase at low T. 
In order to evaluate the glass transition temperature T g , 
Xsg is plotted against N on a log- log scale in Fig. Ub) . 
Strikingly, all the results of XSG give no sign of attaining 
a power law of TV at any finite T we have tried. This 
implies that a true singularity lies only at T = 0. We 
thus conclude that T g = 0; i.e., this system undergoes a 
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FIG. 4: (Color online) (a) Spin-glass susceptibility xsg as 
a function of T. (b) xsg plotted against N^ 1 at T = 
0.100,0.083,0.071, and 0.063, from bottom to top. It does 
not exhibit a power law with respect to N for any finite T. 



zero-temperature glass transition. More detailed finite- 
size analysis remains to be done largely due to lack of 
available system sizes. One interesting point is that a 
hyperbolic lattice does not possess scale invariance as it 
has its own length scale, i.e., the radius of curvature [2rjj ]. 
which may modify the nature of the singularity. 

In summary, we have demonstrated an apparent zero- 
temperature oricntational glass transition in the XY spin 
model on a negatively curved surface. MC simulations re- 
vealed that lowering T makes high-energy bonds gather 
to form clusters, which remain floating with very long 
relaxation times. The long lifetime fluctuation in clus- 
ter distribution follows an Arrhenius-type relaxation at 
low T and the singularity of spin-glass susceptibility 
is expected to arise only at T 9 = 0. These observa- 
tions are consequences of geometry of the surface, where 
curvature-induced frustration in local (and thus global) 
spin configurations yields ground-state degeneracy. 
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